Regularization of the spectral problem for the monolayer graphene with the separable 
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The electronic states in the monolayer graphene with the short-range perturbation asymmetric 
, with respect to the band index are analized. The study was made for the separable in the angular 

momentum representation potential basing on the (2+l)-dimensional Dirac equation. The charac- 
teristic equation for bound and resonance states obtained in the present paper is compared with one 
derived earlier for the same problem with different approach. The momentum representation ap- 
proach used in the present paper allowed us to obtain the satisfactory regularization of the Hadamar 
incorrect boundary problem stemming from the potential singularity. 
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I. INTRODUCTION 



Transport theory in the presence of the resonances and theory of the optical absorption in graphene need nonper- 
turbative analysis of electronic states. This makes it necessary to consider exactly solvable models of defects. One 
' of them is the delta function potential. In the case of the two-band nonrelativistic problem for three-dimensional 
zero- and narrow-gap semiconductors described by the Dirac equation the problem of bound and resonance states 
I . was considered with a use of such potential in the paper The delta function is determined in this model on the 
circumference of the sphere of proper dimension. In the present case it is the circumference of circle. This potential 
has no singularity at r = and is separable in the angular momentum representation. This potential is sepable in 
the angular momentum representation. We take into account possible difference of the perturbation matrix elements 
calculated on wave functions of the upper and lower bands that is equivalent to consideration of both potential and 
mass perturbations. Such asymmetry with respect to the band index can be induced by the local shift of sublattices 
in the vicinity of the point defect. The aim of the present work is a regularization of this incorrect in Hadamar's sence 
problem. 



II. BASIC EQUATIONS 



During the last years much attention was paid to the problem of the electronic spectrum of graphene (see a review 
0). Two-dimensional structure of it and a presence of the cone points in the electronic spectrum make actual a 
comprehensive study of the external fields effect on the spectrum and other characteristics of the electronic states 
described by the Dirac equation in the 2+1 space-time. Short-range potential impurities in graphene were considered 
in works [|| , [3| , [1| . In our works @ , 0] , a new model of the short-range impurities in graphene was considered taking 
into account possible local shift of sublattices. This means that the perturbation must be generically described by 
a Hermitian matrix. We considered it in the diagonal representation. We do not take into account the inter-valley 
transitions. The characteristic equation for the bound states wasderived in [f| within the framework of this model. 
It was understood that singularity of the delta-potential induces a problem of incorrectness of the boundary problem 
that made some regularization necessary. The present work is just dedicated to one possible way of regularization 
of this Hadamar incorrect problem. Notwithstanding a presence of this difficulty, the delta potential is extremely 
popular particularly for use in electronic kinetics and, therefore, deserves thorough analysis. 
The Dirac equation describing electronic states in graphene reads 

(^-ihv F J2 - 7o (m + 5m) vj}j ip = (E -V)ip, (1) 

where vp is the Fermi velocity of the band electrons, 7 M are the Dirac matrices 

7o = o" 3 , 7i = oi, 72 = o- 2 , (2) 
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cr, are the Pauli matrices, 2itivf 2 = E g is the electronic bandgap, ip (r) is the two-component spinor: 

^( r )=f{£iY (3) 



9 W 

. The electronic gap can appear in the graphene monatomic film lying on the substrate because of the sublattices 
mutual shift Q. The spinor structure takes into account the two-sublattice structure of graphene. 5m (r) and V(r) 
are the local perturbations of the mass (gap) and the chemical potential. A local mass perturbation can be induced 
by defects in a graphene film or in the substrate [|[ . We consider here the delta function model of the perturbation: 

5m (r) = —b5{r — r ), V{v) = —aS(r — r ), (4) 
that can be re- written in the form of the band- asymmetric potential 

Vx(r) =V?5{r - r ), V 2 (r) =V 2 °S(r - r ). (5) 
The parameters V°, V 2 ° , a and b are related as follows 

V? = - (a + 6) , V 2 ° = - (a - b) . (6) 

Here r and ro are respectively the polar coordinate radius and the perturbation radius. Such short-range perturbation 
was used in the (3-(-l)-Dirac problem for narrow-gap and zero-gap semiconductors in The perturbation forms the 
diagonal matrix 

diag{V?X) ( 7 ) 

The delta function perturbation is the simplest solvable short-range model. Finite radius ro plays a role of the 
regulator and is necessary in order to exclude deep states of the atomic energy scale. The finite perturbation radius 
ro leads to the quasi-momentum space form-factor proportional to the Bessel function that justifies our neglect of 
transitions between the Brillouin band points K and K' [l| . 

Let us introduce the two-dimensional Fourier representation of the two-component wave function ([3l 

with the Fourier-components / p and g p determined by the inverse Fourier transforms: 

/ p = J dxdye- ikr f (r) , g p = J dxdye^g (r) (9) 
The inverse transform can be written as a combination of the Hankel transform and the angular Fourier series 

oo oo 

U = f(p,e)= Y, J?v " e /n(p), gp = g(p,o)= ]T i n e in9 g n (p) (10) 



fn (p) = / drrj n (r) J„ (pr) , g n (p) = / drrg n (r) J n (pr) (11) 
Jo Jo 

where /„ (r) is the angular Fourier component: 

/(r,0) = -= £ f n (r)e me , g(r,9) = -= £ g n (r) e™ e (12) 



'2n ^ V2tt 

n — — oo v n— — oo 
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Making the Fourier transform of the Dirac equation ([T]), we consider preliminary two terms of it, which can be 
considered as rather complicated. The kinetic term reads 

-i(^ + ^)^(r)=(| : +|)^ ( r r )). (13) 

Therefore, we can write for the upper and lower components of the kinetic term in the momentum representation 

fax + ip y ) g (p) \ _ ( pe t9 g(p,9) \ , . 

(p x -ip y )f(p) J ~ \ P e- l9 f{v,e) 



Substituting (fit)]) into (TUT) , we obtain: 



-iJZZ-oo^fn-lip) 



(15) 



The potential Fourier transform Vi (p) can be expanded into a series in the case of the circular symmetry: 

oo 

V i (\p-p'\)=-K e me -V?(p,p>), (16) 



V?(p,p')= / dr-rJ n {pr)Vi(r)J n {p'r) (17) 
Jo 

where J n (pr) is the Bessel function, i = 1, 2. Thus we obtain the integral equations: 



(E - m) fj (p) + ip 9j (p) - j dp'- p'fj (p') V( (p',p) = 0, (18) 

(E + m) 9j (p) - ipfj (p) - J dp'- p' 9j {p') Vi (p',p) = 0, (19) 

where we have put n = j — 1/2, where j is the pseudospin quantumnumbcr, j = ±1/2, ±3/2 .... In opposite to 
the relativistic theory, this quantum number has nothing to do with the real spin and indicates a degeneracy in the 
biconic Dirac point. These equations have a symmetry: 

ft O g h E -> -E, j - 1/2 -> j + 1/2, a -► -a. (20) 

We have introduced the notations: 

/n=j-l/2 = fj: 9n=j+l/2 = 9j (21) 

Zero-radius [lOj and separable potentials [Tl| are popular in the nonrelativistic scattering theory. However, the Dirac 



equation is extremely sensitive to a singularity of the potential 12j . The singularity of the delta function potential 
can be regulated putting the delta function to the circular support (|4|) [l| 

Substituting ([7]) into (|16p we obtain the separable in the angular momentum-momentum modulus representation 
potential: 

Vi(P,P')=4(p)4(p'), (22) 

where 



4 (P) = yfioV? ^--1/2 (^o), «2 (P) = Jn>v]jj+i/2(pro),i = 1, 2 (23) 
Equations (|T8|) . (fl9|) become degenerate and can be written as follows: 



oo 



(E - m) fj (p) + ip 9j (p) - v{ (p) / dp'p'ft (j>') «J (p') = 0, (24) 

Jo 

/•oo 

(£ + m) 5j (p) - tp/j (p) - ^ (P) / Wft (p') «| (p') = 0. (25) 



o 



III. CHARACTERISTIC EQUATION 



Introducing the functions 

Fj (E) = / dppfj (p) v{ (p) , Gj (E) = I dppg 3 (p) v\ (p) , R(p) = (p 2 + m 2 - E 2 ) _1 , 



jo Jo 
we obtain the homogeneous algebraic equations set for F and G: 
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Fj I 1 



+ (E + m) j dppR (p) (v{ (p)) J - iGj j d PP 2 R (p) ^ (p) (p) = 0, 
jf dpp 2 i? (p) wj (p) ^ (p) + \l + {E-m)j dppR ( P ) (v> 2 (p)) ^ = 0. 



The solvability condition for this equations set gives the characteristic equation: 

2' 

) 
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1 + (m + E) I dppR (p) (v{ (p)j • 1 + (m - E) f dppR (p) ^ (p)^ 



n 2 



dpp i? (p) (p) ?4 (p) 



Using the well known formula [1 

dx 



X 2 +Z 2 



J M (6a:) J, (ex) = (-1)" z^ v+2n I v (x) (x) 



we write it in the form: 



[1 + (to + E) roV°I^ 1/2 (wo) i^-1/2 [1 + (m - E) V 2 °I j+1/2 (nr ) K j+l/2 (nr )] 

= (m + E)(m- E)V?V 2 °r 2 I*_ 1/2 ( Kr °) K j+i/2 ( Kr °) > 

where k 2 = (m 2 — E 2 ) , /„ (x) , K n (x) are the modified Bessel functions. Making use of the identity [IH 

I v (x) K v+ i (x) + I v+ x (x) K u (x) = 1/x, 
we obtain the characteristic equation: 

« [ij-i/a ( Kr o) («ro) + ^+i/2 («ro) ^-1/2 («ro)] +(m + E) V°Ij-yi («x ) ^-1/2 («ro) + 

(to - £?) V A 2 / J+ i/2 (w ) (wo) = 0. 

Using the relations ([6]) , we can write the characteristic equation in the form: 

K [lj-1/2 ( Kr o) Kj+1/2 ( Kr o) + Kj-1/2 ( Kr o) ( Kr o)] = 



[(to - E)(a- b) I J+ i/2 (nr ) K J+1/2 (/cr ) + (a + b) (to + E)I^ x /2 («ro) Kj-1/2 ( Kr o)] ■ 
This equation has to be compared with one derived using the different approach in 0] : 



P (jj-1/2 {pro) H f+i/2 (P r o) - J3+1/2 {pro) H f\/ 2 (P r o)) 



tan (V« 2 - b 2 ) 
Va 2 - b 2 



E — TO 

E + TO 



(a - b) J j+1 / 2 {pro) H j+i/2 (P r o) + 



jjrz^ ( a + b ) {Pro) Hf} 1/2 (no) 
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where (z) is Hankel's function, p = \/E 2 — m 2 . Making the analytic continuation from the case of the band states 
to the opposite one of bound states E 2 < m 2 we obtain the equation 

K [ I j-l/2 («fo) K J+1/2 («r Q ) + ^-1/2 («m) («r )] = 



tan (Va 2 - b 2 ) 



[(m — E){a— b) (kxq) #i+i/2 (wo) + (a + 6) (m + E)I J _ l/2 (nr ) K j _ 1 / 2 (nro)] • (35) 



Va 2 - 6 2 

We see that the only distinction between the formulae (|33[) and (|35|) is a presence of the factor 



, s tan(Va 2 -6 2 ) , , 

T (a, b) = \ ' 36 

V ^ Va 2 - b 2 V ' 

in the equation derived in [7J . They are identical in the limit 

a 2 -b 2 — >0, T{a,b) — > 1. (37) 

This limit can be reached near the lines a 2 — b 2 = 0, or when a and b are small. The reason of this distinction is 
the following. The delta function potential in the Dirac equation creates a boundary problem, which is incorrect 
in the Hadamar sense [l5j]. Such a problem needs a regularization. Distinct regularization procedures were used 
in the present work and in Q. Singular feature of the potential was preserved in Q, a regularization was carried 
out due to very special join of partial solutions, while in the present work, a regularization took place due to the 
Fourier transform smoothing property. In particular, the potential takes the separable form with the nonlocal kernel 
v i (p) = \/ r Jj+1/2 {P r o ) (see eq. ([23")) h Just this nonlocality plays the role of regularizator. Equations (1531 and 
([3"5j) have asymptotically similar solutions , when (|3"T1) is valid. Their sets of solutions are quite different outside this 
region of parameters a and b. Moreover, equation ([35]) has significantly richer set of solutions, than equation ([33]) 0] 
due to periodicity of tangent present in (1551) . We come to conclusion that the rich set of solutions given by ([55)1 is 
the artifact of the zero-support interaction. Besides, the annular well model analyzed in [16| does not reproduce all 
of this rich set of solutions even in the limit of zero radius. Therefore, the equation obtained in the present work can 
be considered as a correctly regularized one. 

IV. CONCLUSION 

Non-relativistic problem of the electronic spectrum described by the Dirac equation is considered in the case of the 
two-component short-range perturbation potential with a help of the Fourier transforming. This approach is shown 
to make regularized the incorrect in the Hadamar sense problem. 
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